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Abstract: After the prediction that strong enough optical activity 
may result in negative refraction and negative reflection, more and more 
artificial chiral metamaterials were designed and fabricated at difference 
frequency ranges from microwaves to optical waves. Therefore, a simple 
and robust method to retrieve the effective constitutive parameters for chiral 
metamaterials is urgently needed. Here, we analyze the wave propagation in 
chiral metamaterials and follow the regular retrieval procedure for ordinary 
metamaterials and apply it in chiral metamaterial slabs. Then based on 
the transfer matrix technique, the parameter retrieval is extended to treat 
samples with not only the substrate but also the top layers. After the 
parameter retrieval procedure, we take two examples to check our method 
and study how the substrate influences on the thin chiral metamaterials 
slabs. We find that the substrate may cause the homogeneous slab to be 
inhomogeneous, i.e. the reflections in forward and backward directions 
are different. However, the chiral metamaterial where the resonance ele- 
ment is embedded far away from the substrate is insensitive to the substrate. 
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1. Introduction 

Chiral metamaterials (CMMs), which are artificial materials that lack any planes of mirror sym- 
metry, possess strong ability to rotate the plane of polarization of electromagnetic waves. This 
ability is called optical activity which is of great interest to many areas of science, for exam- 
ple, analytical chemistry and molecular biology. Recently, it was predicted that strong enough 
optical activity may result in negative refraction Q]-[6] and negative reflection (7). After these 
predictions 0]-[7), a lot of artificial CMMs were designed and fabricated at different frequency 
ranges from microwaves to optical frequencies [ 8-22 1 . The basic physical quantities to char- 
acterize the optical properties of CMMs are the effective constitutive parameters. Therefore, a 
simple and robust method to retrieve the effective parameters for CMMs is urgently needed. 

Parameter retrieval I23H281 is a basic and important technique to obtain the electromag- 
netic properties of the effective media. The effective media are usually considered to be ho- 
mogeneous when the size of structure is much smaller than the wavelength A. The retrieval 
parameters, e and /i or n and Z, are well defined and can be determined from reflection and 
transmission coefficients (S parameters) (2314281 . We follow the similar procedure (23l to do 
the parameter retrieval for CMMs as it was done in ordinary metamaterials. The procedure of 
the retrieval for free-standing CMM slabs was briefly shown in our previous paper (29). How- 
ever, in most of the experiments 1 17-22], the CMMs are fabricated on a substrate. Therefore, a 
retrieval procedure for samples with substrates should be developed. This is an important devel- 
opment for obtaining the electromagnetic properties of CMMs from experimental data. Kwon 
et al. Il30l presented a retrieval procedure for CMMs with a semi-infinite substrate. However, 
experimentally the thickness of the substrate is finite. Here, we present a simple and robust 
retrieval method based on the transfer matrix technique. It can deal with samples not only with 
substrates but also with top layers. 

After the parameter retrieval, we give two chiral metamaterial examples: 1) twisted-crosses 
and 2) four-folded rotated i2-particles. For the twisted-crosses CMMs, the resonance elements 
are very close to the substrate. The existence of the substrate, therefore, has a huge influence on 
the thin CMM slab. It causes the resonance shifts to low frequencies and induces the homoge- 
neous slab to be inhomogeneous. In order to correct the inhomogeneity, we need to introduce a 
top layer, the same material as the substrate, to cover the CMM slab. The top layer will red-shift 
the resonance further. For the four-folded rotated £2-particles CMMs, the resonance elements 
are embedded far away from the substrate. The influence of the substrate can be neglected. In 
the meanwhile, profiting from the retrieval method and the precise simulation, we fit the re- 
trieval results using the analytical formulas deriving from the ii-particle resonator model. They 
agree with each other very well. 

The following sections are arranged as follows. In Sec.|2] we first obtain explicit analytic ex- 
pressions of the retrieval parameters for free standing CMM slabs, then we extend this method 
to deal with the samples with the substrate and top layers. In Sec. [3] we use the parameter re- 
trieval methods to study two CMM examples and obtain the influence of substrate/top layer on 
thin layer CMM slabs. Finally in Sec. [4] we present our conclusions. 



2. Retrieval method 



2.1. Eigen waves in chiral materials 

For reciprocal (Pasteur) chiral materials |3 111321 , the strength of the cross coupling between the 
magnetic and electric field is characterized by the parameter K called chirality. The constitutive 
relations in chiral medium are usually written as 11331 : 

£ £ iK/co \ /E\ 

-iK/c jUoM / V H / ' 

where £q and /Iq are the permittivity and permeability of vacuum, £ and ji are the relative 
permittivity and permeability of the medium respectively, cq is the speed of light in vacuum. 
Assuming all fields to be plane waves, the source free Maxwell equations take the form: 

kxH= coD, kxE = oB, (2) 

where k and co are the wavevector and the angular frequency of the plane wave in the chiral 
medium. Inserting (Q~|) into (0, we can get the eigenfunction of the electric field E propagating 
in the chiral medium: 

kx (kxE) = -^(£ J Li-fC 2 )E-2/K-A:o(kxE), (3) 

where ko = (0 /co is the free space wavevector. For simplicity but without loss of generality we 
can assume k = kz and obtain the eigenvectors and eigenvalues as: 

E ± (z) = (x±iy)E Q e ikz , (4) 

k ± =k Q (n±K), (5) 

where n = y/ejl. There are two eigenmodes in chiral medium. One is the right circular polar- 
ization (RCP/+) wave and the other is the left circular polarization (LCP/-) wave. Any wave 
propagating though the medium can be uniquely decomposed into RCP and LCP and there is 
no conversion between RCP and LCP as they propagate in the chiral medium. Then we can 
define the index of refraction for RCP/LCP waves as: 




n± =n±K. (6) 

If K > n, n_ will be negative. This constitutes an alternative approach to realize a negative 
refractive index as proposed earlier by Pendry (TJ. 

A chiral medium has two important properties. One is called optical activity which charac- 
terizes the rotation of the polarization plane of a linearly polarized light as it passes through a 
chiral medium. Mathematically, it's defined as the polarization azimuth rotation angle of ellip- 
tically polarized light: 

= i[arg(r + )-arg(7L)], (7) 

where T + and T are the transmission coefficients for RCP and LCP waves. The other prop- 
erty is circular dichromism which arises from the different absorption for RCP and LCP. It 
characterizes the difference between the transmissions of two polarizations: 

n _i -i r ir+l 2 -|r-r\ (S) 



Artificial CMMs with large 9 and small 77 are very important for applications. 
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Fig. 1. Schematics of the transmission and reflection coefficient of a standalone CMM slab. 



2.2. Parameter retrieval for CMM slabs without substrate 

A slab of reciprocal (Pasteur) chiral material looks exactly identical for waves propagating 
in forward or backward direction. However, for a normal incident circular polarized wave 
E± (z) = (x ± iy)Eoe' kz along the z direction the reflected wave will be E' ± (z) — (x ± iy)EQ e~' kz 
along the — z direction. The y component is always 90° advanced(+) / retarded(-) with respect 
to the x component when the wave vector changes to the opposite direction. Therefore, the 
polarization of the reflected wave is reversed based on our current definition of handiness in 
relation to the wave vector, i.e. the reflected wave of RCP (LCP) wave is LCP (RCP). Now we 
consider a circularly polarized plane wave normally incidents onto a CMM slab with thickness 
d, refractive index n±, and impedance Z = \/jJ./e, as shown in Fig.Q](note that the impedance, 
Z's, of RCP and LCP waves are identical |3T1 ). The amplitude of the incident wave is assumed 
to be unity; the amplitudes of the transmitted and reflected waves are T± and R+ respectively; 
those of the forward and backward propagating waves inside the slab are T' ± and R^. Applying 
the condition of continuity of tangential electric and magnetic fields at z = and z = d,we have 
the following equations: 



l-R T = ± - i (9b) 



l+R T = Ti+R 1 ^, (9a) 
Tl - R'^ 
~Z 

T' ± e ik±d + R'^e- ik T d = 7± , (9c) 
jf e ik±d _ R i -ik T d 

— ^ = T± , (9d) 

where k± is the wave vector as defined by Eq. (0. Then eliminating T' and R' in Eqs. (O, we 
obtain the transmission and reflection coefficients: 

4Ze' nk ° d e ±iKk ° d 

T± = (l+Z) 2 -(l-Z) 2 e 2 »'M' (10a) 



(l-Z 2 )^ 2 '"^-!) 
(l+Z) 2 -(l-Z) 2 e 2 '"M ■ 



Eq. (fTOb) shows that R + = R . This can be understood intuitively as follows: We know that 
the polarizations of RCP and LCP are reversed after being reflected. For the forward incident 
wave, the optical path of RCP (LCP) is n+d (ri-d), where d is the distance that the wave has 
traveled; For the backward reflected wave, the optical path is n-d {n+d). Therefore, the total 



optical paths for RCP and LCP are all n + d + ri-d. They pass through the same optical path 
for the reflected waves, so the loss should be the same for both of them. In addition, the same 
impedance of the RCP and the LCP gives the same reflections for both circular polarizations 
on the surface. Therefore, both the reflections of the single slab should be the same. 

Inverting the above Eqs. ( [Tol l, the impedance Z and refractive index n± are given below: 



Z=± f(i-Ry-T + T_ > (lla) 

^ = ^{H^d-^m±2mn}, (lib) 

where m is an integer determined by the branches. The correct sign of the square root in Eq. 
( II 1 al l and the correct branch of the logarithm in Eq. dl lbb are selected according to the follow- 
ing conditions: 

Re(Z) >0,lm(«) > 0, (12) 



which are required by energy conservation and causality 1231 . Then the remaining retrieval pa- 
rameters can be obtained with the following relations: n = (n + +«_)/2, K = (n + — n_)/2, £ = 
n /Z, and = nZ. 

Actually, we have another route to get the retrieval results. If we redefine T = y/T + T- and 
R = R^, according to Eqs. (flOt . T and R can be written as: 

4Ze ink ° d 

T= (l+Z) 2 -(l-Z) 2 eW (13a) 

(1 -Z 2 )(e 2i " k ° d - 1) 
R= (l+Z) 2 -(l-Z) 2 eW (13b) 

These expressions are exactly the same as those for the ordinary nonchiral materials. There- 
fore, we can use the traditional retrieval procedure for the nonchiral metamaterials to get the 
refraction index «, the impedance Z, the permeability e and the permittivity jj, 1231 . Then from 
Eqs. (fTOb), we can obtain 

K=——\nC ± ) = —H + . ), (14) 
2k d y T-' 2k d 

where <p± are the phases of T±. Then we immediately get the refraction indexes n± = n±K. 
From Eqs. (fl4l . we have 

<b+ — (b-+ 2mK 

Re ^= 2M ' (15a) 
= m\T-\-m\T + \ 

2k(jd 

where the integer m is determined by the condition of — n < <j> + — 0_ +2mn < % for one unit 
cell. We note that the real and imaginary parts of K relate to the azimuth rotation angle 9 and 
the circular dichromism rj respectively. 

2.3. Parameter retrieval for CMM slabs with substrate 

The retrieval procedure discussed above identifies the scattering amplitudes of the CMM slab 
with those of an equivalent homogeneous slab. The retrieved effective parameters of the CMM 



are the parameters of this equivalent homogeneous slab and given as functions of the scattering 
matrix of the CMM. In this form, the retrieval procedure is only applicable if the CMM has the 
same symmetry as the equivalent homogeneous slab. For instance, the reflection amplitudes for 
a homogeneous chiral (Pasteur) slab are identical for incidence from the front or back interfaces. 
For a CMM on top of a substrate (for an example see Fig 2) the scattering amplitudes do 
not necessarily have this property. In this case, we cannot directly apply the above retrieval 
procedure to the measured or simulated scattering amplitudes of the sample to obtain effective 
parameters. 

In the experiments and simulations (for instance, the commercial software CST Microwave 
Studio ll34l ). we can easily obtain the total scattering matrix S tota i of the CMM/substrate slab 
which relates the incoming waves to the outgoing (scattered) waves. It is well known l35ll that 
any scattering matrix S corresponds to a transfer matrix M which relates the amplitudes of the 
in- and outgoing waves in front of the slab to those behind the slab: 

c (r t'\ (t-rr'/t' r> lt'\ 

S =[t ,'J ; M =( -r/t> l/t>)> (16) 

where t and r are the transmission and reflection amplitudes; the quantities without (with) prime 
(/) refer to wave propagation in forward (backward) direction. 

Now we consider a structure of substrate/CMM/top layer as an example and assume the pos- 
itive (forward) direction to be from the top layer to substrate layer. Measurement or simulation 
afford the total scattering matrix S tota i of the stack in terms of total transmission and reflection 
amplitudes. Using ( fTol l we obtain the corresponding total transfer matrix M tota i, which factor- 
izes into a product of the transfer matrices of the individual layers [35): 

M tota i = M su bstrateMcMMM t0 player ■ (17) 

The homogeneous substrate and top layers are assumed to be known. Their transfer matrices, 
M su bstrate an d M t0 pi ayer , respectively, are just the well known transfer matix of a homogeneous 
slab and can be calculated analytically (as the product of the transfer matix across the interface 
from vacuum into the homogenous slab, the diagonal transfer matrix (phase & attennuation) for 
the propagation within the homogeneous slab, and the inverse of the interface transfer matrix 
again to step out of the slab back into vacuum): 



M 



z+i z-i 
zfi ih 

2Z 2Z 



Jnk d Q 
e -i"k>d 



1+Z 1-Z 

2 2 

1-Z 1+Z 

2 2 



(18) 



where the impedance of Z, refractive index of «, and the thickness of d are all known parameters 
of the substrate or the top layer. Note that the transfer matrix of each layer is that of the slab 
in free space, i.e. the interface of two media can be considered as two media separated by a 
vacuum slab with zero thickness. 

Solving the above equation $1% , we can get the effective transfer matrix of the CMM slab 
standing alone in free space: 

McMM = M s ~ubstrate M totalM to p layer . (19) 

Having extracted the effective transfer matrix of the CMM slab alone, we can use the relation 
( [Tol l again to obtain the S parameters of the free standing CMM slab. Then the retrival procedure 
given in Sec. |2.2| can be applied to the isolated S parameters of the CMM slab. And following 
this way, we can obtain the retrival parameters of CMMs with substrate and/or the top layer. 



3. Examples of the CMM slabs with and without substrate 



In this section, we will give two examples to demonstrate that a physically meaningful retrieval 
can be obtained from scattering amplitudes measured for the combined system of CMM and 
substrate/top layer and to show how the substrate influences the retrieval parameters of the 
CMM slab. For this purpose, we choose the twisted-crosses CMMs which had been fabricated 
in the optical region by Decker et al. ll2D and the four-folded rotated H-particle CMMs which 
is closest to the standard chiral model. 



3.1. Retrieval results of the twisted-crosses CMM slabs with and without substrate 

air/CMM/air air/top layer/CMM/substrate/air air/CMM/substrate/air 
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Fig. 2. (a)-(c) are the Schematics of the twisted-crosses CMM slab, (a) free standing, (b) 
with substrate and top layer, and (c) with substrate only. The amplititudes, (d)-(f), and 
the phases, (g)-(i), are also shown correspondingly. The prime (/) denotes that the wave is 
incident from the opposite direction. 

The structure of the twisted-crosses CMM is shown in Fig. [2] It's composed of two gold 
crosses rotated against each other by an angle of 22.5°. Each arm of the crosses is a short 
rectangular wire with the dimension of 25 nm x 56 nm x 315nm. The separation between the 
two crosses is 37 .5 nm. The gold structure is embedded in a lossless dielectric with the refractive 
index of n = 1.41 as shown in Fig. |2|a). The structure is fabricated on the glass substrate 
(n=1.45) with the thickness of 500nm as illustrated in Fig. [2c); Optionally a thin glass layer 
(n=1.45) with the thickness of 200 nm can be placed on top of the embedded gold structure, 

CO 2 

see Fig.|2b). The gold is modeled by a Drude model permittivity e(fi)) = £» — m ( m + jm ) with 
plasma frequency co p = 2n x 2159THz, collision frequency co c = 2n x 25 THz, and £„ = 9.07. 

The lower two rows in Fig. [2] show the circular polarization scattering amplitudes (magni- 
tudes and phases) of the CMM slab obtained from simulation. For the free standing slab (left 
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Fig. 3. The retrieval parameters £, fl, K, n + , tj, and 6 for CMM slabs. (Left) free- 
standing, (Middle) with substrate and top layer, and (Right) with substrate only. These re- 
sults were extracted using 



column) the simulation directly yields the scattering amplitudes of the CMM slab; for the CMM 
on substrate (right) and CMM on substrat with top layer (middle column) the shown CMM scat- 
tering amplitudes were extracted using (TT9b derived above. T (R+-) and 7++ (R |_) are the 

transmission (reflection) of LCP and RCP, respectively. (In order to specify the polarization 
conversion, for T a i, and R a b, b refers to the incident wave and a to the transmitted/reflected 
wave.) As can be seen in Figs. [2] the reflection amplitudes for RCP and LCP are exactly the 
same, which confirms the previous analysis. 

The scattering amplitudes of the homogeneous chiral slab are the same for propagation 
though the slab in forward or backward direction. The simulated scattering amplitudes for the 
free standing CMM shown in Figs.|2ja,d,g) also have this property. We can directly apply the 
retrieval procedure described in section |Z21 and obtain the effective medium parameters. The 
results are given in Figs. [3] 

For the CMM on susbtrate no meaningful retrieval can be obtained by applying the retrieval 
procedure described in section l2~2l directly to the total scattering amplitudes; The presence of 
the substrate leads to additional scattering not originating from the CMM. It also makes the 
reflections in forward and backward direction unequal - an asymmetry which can never be 
reproduced by the homogeneous chiral slab. Only after extracting the scattering amplitudes of 
the CMM slab itself from the simualtion results using Eq. (TT~9b derived in section [231 we can 
apply the standard retrieval procedure described in section l2~2l to the extracted CMM scattering 
amplitudes and obtain reasonable effective medium parameters for the CMM slab. 

Even after extracting the CMM scattering amplitudes we still observe some weakly asym- 
metric reflections R + ^ R' + _ in forward and backward direction for the CMM on substrate in 
Figs.|2f,i). The extraction derived in section 1231 does only deal with the propagating modes 
in the system. However, there is also interaction between the CMM slabs surface regions with 
their environment mediated by evanescent modes. This interaction makes the surfaces of the 
CMM slab on substrate unequal and prevents its approximation by an equivalent homogene- 
nous effective medium slab. 

There are two possible solutions to this problem: First, we could follow the method intro- 
duced in Ref. 22 and approximate an isotropic reflection R av as the geometric average, 

R av = \/RR^, (20) 

where, R (/?') is the reflection when the incident wave is from the top/vacuum (bottom/substrate) 
side of the CMM slab. This method has been used successfully if the anisotropy was small 
enough 11271 . (Note that the inhomogenization is not inherent property of the twisted-crosses 
CMM slab. It's only the boundary effect of the thin layer sample. If a sufficiently transparent 
sample goes to bulk material with infinity thickness, this boundary effect becomes negligible.) 
Using this method we obtain the retrieval of the single sided CMM on sustrate shown in the 
right column in Fig. [3] 

Second, we could cover the top side of the CMM with a second layer of substrate. This layer 
does not need to be very thick as the interaction via the evanescent modes with the surface 
layers of the CMM decays quickly. As can be seen in Figs.|2e,h) this approach lead to perfectly 
symmetric reflection amplitudes and allows to apply the standard retrieval procedure described 
in section l2~2l to the extracted CMM scattering amplitudes to obtain the effective parametes of 
the CMM slab. The results for this case are show in the center column in Fig. [3] 

Comparing the results with and without substrate, we can see that the substrate has two 
main influences on the CMMs. One is lowering the frequency of the resonance because the 
substrate increases the capacity C of the resonator and (Or ~ 1 / VLC. The other is the substrate 
may induce the homogeneous CMM to be inhomogeneous. For the structure with substrate 
only, the two crosses are in the different environments. For the cross near the substrate, the 



evanescent field leak to the substrate while for the cross near the air, the field leaks to the 
air, which makes the two substantially identical crosses unequal physically. Even the substrate 
contribution has been deducted off mathematically by using Eq. dT9l ). the CMM slab still shows 
the inhomogeneous property (e.q. ^ /?+_)■ 

3.2. Retrieval results of the four-folded rotated Q.-particle CMM slabs with and without sub- 
strate 




Fig. 4. Schmetics of the four-folded rotated fi-particle CMM. The inner and outer radius 
of the loop are r = 5 flm and R = 5.5 jlm. The length of each straight wire is L = 8.08 flm. 
d = 0.5 /im. a = 32 jlm. b = Hjxm. Each fl-particle locates at the center of each quarter 
of the unit cell. The silver fl-particle is embeded in the polyimide with e = 2.5 and loss 
tangent S = 0.03. The silver is charactered by the Drude model with the surface plasmon 
(Op = 13.66 x 10 15 rad/s and the collision frequency co c = 2.73 x 10 13 rad/s. 

Not all the chiral metamaterials are so sensitive to the influence of the substrate. In the fol- 
lowing, we will give such an example — four-folded rotated £2-particle chiral metamaterials. 
The structure and the dimensional parameters are shown in Fig. [4] The incident wave is per- 
pendicular to the CMM slab. 

The amplitude and phase of the reflection and transmission are shown in Figs. |3a-d). The 
left is with 50jtim substrate (n=1.27) only and the right is for the free standing sample. There 
is a strong resonance at around 1.87THz where both the amplitude and phase of the transmis- 
sions of RCP and LCP are different, which indicates very strong optical activity and circular 
dichromism. The interesting is that for the structure with substrate only, the extracted reflec- 
tions from different directions almost overlap with each other, e.g. R = R', even they have a 
large difference before the extraction using Eq. (19[ . Comparing the results between with (Fig. 
|5Ja,c)) and without (Fig.|5Jb,d)) substrate, we see that the substrate almost has no influence on 
the £2-particle CMMs. That is because the resonance elements, the gap of the circle magnetic 
resonator and the short wires for the electric resonance, are located in the middle of the unit 
cell. They are far away from the unit cell surface, in other words, they are far away from the 
substrate. The evanescent field leaked into the substrate is very little. Therefore, this kind CMM 
is not sensitive to the existence of the substrate. 

Figs. HJe-j) are the retrieval results of the optical parameters £, ji, and K. All these curves 
behave well. We fit these curves using the analytical formulas given by Eqs. (|29l derived from 
the H-particle resonator model (see Appendix [A}. For with substrate only, the fitted param- 
eters are as following: e b = 3.1738, fi b = 0.9799, coq = 1.8651THz, y = 0.05519a)o, £2 £ = 
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Fig. 5. (a)-(d) are the simulation results of the amplitude and the phase of the reflection and 
transmission for the four-folded rotated fl-particle CMM slabs. They are extracted using 
1191 . R 1 is the reflection for the incident wave from the opposite direction, (e)-(j) show the 
optical parameters e, jl, and K. Re/Im denotes the real/imaginary part of the value. The 
subscripts s and f denotes the retrieval results from the numerical simulation and the fitting 
results using Eqs. J29b respectively. The left is with substrate (n=1.27) only and the right is 
for the free standing sample. 

0.1537, £2^ = 0.0627, andQ. K — 0.0986. For without substrate, they are very close to those with 
substrate. e,, = 3.1736, ^i = 0.9798, fflo = 1.8713THz, 7= 0.05463 too, £2 £ =0.1560, £2^ = 
0.0625, and D. K = 0.0993. From these parameters, we can again see that the substrate has very 
little influence on the four-folded rotated £2-particle CMMs. The substrate only lowers the fre- 
quency of the resonance by 0.3%. 



In Figs. |3e-j), we can also see that the retrieval results from numerical simulations agree 
very well with the analytical results given by Eqs. (1291 ). It verifies that the ft) dependences of the 
optical parameters are valid in our previous work about realizing the repulsive Casimir force 
ll36l using chiral metamaterials. Note that from Eq. ( 1251 ). the strength of electric polarization 
and magnetic polarization can be modified independently by the length of the straight wires and 
the area of the circular loop respectively. £2 £ > Q.^ here is just one case. But no matter how to 
change the dimensional parameters, the relation Eq. ( f3Qb roughly always holds in such designs 
composited by the passive material only. 

4. Conclusion 

In summary, we have analyzed the wave propagation in chiral metamaterials and followed the 
similar procedure that was done in ordinary metamaterials to do the parameter retrieval for 
CMM slabs. Based on the transfer matrix technique, we extended the parameter retrieval tech- 
nique to be able to treat samples with substrates and extra top layers. We studied the influence 
of the substrate and top layers on the thin chiral metamaterial slabs. We found that the substrate 
could lower the frequency of the resonance and induce the homogeneous chiral metamateri- 
als to be inhomogeneous. The sensitivity of the influence relates to the strength of the leaked 
evanescent field in the substrate. We also fitted the retrieval results using analytical expressions 
derived from H-particle chiral model based on the effective LC circuit mode. We found they 
agree with each other very well. 
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A. Derive the forms of the constitutive parameters from the i2-particle resonator model 



Fig. [6] gives the schematics of the £2-particle resonator. The structure consists of an open 
circular loop and two short wires. The area of the loop is S = nr 2 and the length of each wire 
is I. The wires are perpendicular to the loop and connected to the ends of the open loop. If the 
i2-particle is in a homogeneous external field, the driving electric potential can be written as: 




Fig. 6. Schematics of the single right-handed (left) and left-handed (right) fl-particle res- 
onators. 



U = 21E ±LK}SH , 



(21) 



where ± signs correspond to right-handed and left-handed helix resonators, The hat-dot denotes 



the first time derivative. Applying the effective RLC circuit model, we have 

Li+^+RI = U, I = q. (22) 

Let 0£ = j, /? = jUo j> 7= f, and co^ = Then assuming the external fields are plane waves, 
~ e~' mt , and considering stationary solutions, we get the following equation: 

(-co 2 -i(Qy+(Qo)q = aEo±icoBH , (23) 

then, we have 

a icoB 

q = r~. -^E ± , . , 2 H . (24) 

-co 2 -my+co^ -co- -ig>y+(Oq 

Therefore, the electric dipole, p = q\, and the magnetic dipole, m = ±/A = ±qA, can be written 

as: 

P = a a \ . E Q + /'f H , (25a) 

= TfmaA g + fl^pA g (25b) 

where, aA = 1 = lh, and A = AIq; lo is the unit vector along the direction of the wires, i.e., 

the up direction. The electric polarization and the magnetic polarization are defined as P = ^£ 

and M = respectively, where the summation is done in the unit cell and Vo denotes the vol- 
ume of the unit cell. Assuming there are N resonators in one unit cell, the electric polarization 
and the magnetic polarization can be written as: 



P =f73-^E+^ 7 ^-H, (26a) 



N al £ I N ±iiopl 



M =f T'fflttA E+ N mp A u (26b) 

v o az—ar—iwy Vo <ak— eo z — imy 

Here, the directions of 1 and A are merged into Eq and Hq to form the vectors of E and H. 
Inserting Eqs. d26b to the relation D = eqE + P and B = jUoH + fioM, we have 

D = £0 E + , gf y». E + f /V ° H, (27a) 

B = MoH+ T f W( ; AjV/y0 E+ H. (27b) 

Comparing Eqs. < f27b with Eq. ((TJ, we can obtain the relative permittivity e, permeability jU, 
and the chirality k: 

8 = 1 + , aW/ 2 V ° £0 , (28a) 
(Oq - co 2 - my 

co 2 BAN/V Q 

£i=l + — = r ; , (28b) 

cofi - co 2 - my 

_ ±coBlc N/V _ ±(oii c aAN /V ... , 

K— — 5 5 — : — — — o t — : • (Z5CJ 

coQ-co 2 -my (OQ-(D--my 

Eqs. d28t are valid for the case that the metal chiral structure stands in free space. If it's merged 
into a background material with and /if,, Eqs. (f28l > should be modified as follows: 

e = ei , + ^ / . , (29a) 

co ( j - co 2 - my 



Q. a co 2 

M = M;,+ 2 2 ■ ■ (29b) 
at 1 -my 

where £/, is usually larger than one; /I/, is usually very close to one. D. £ , D.^, and Or- are the 
coefficients of the resonance terms in e, u, and K, i.e. i2 £ = , , alN ■> , i2 u = and = 

= ^y"^^ i where co is the speed of light in vacuum. They describe the strength of the 
resonance. The expressions of the frequency dependence of e and jj. given by Eqs. (|29al > and 



d29b| ) are the same as those in traditional metamaterials 13711 . The linear co dependence of K is 
the same as Condon model for the homogeneous chiral molecular media l38l . This is a general 
feature of natural optically active materials l39l . 

Because a = I and j3 = jtiof - we then can obtain the following relation: 

i4 = £2 £ £2 M , (30) 



which seems to limit the possibility of obtaining large Q K in these designs composited by the 
passive materials only. 



